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Abstract 

Background: Although many molecular regulators of adipogenesis have been identified a comprehensive 
catalogue of components is still missing. Recent studies showed that the retinoblastoma protein (pRb) was 
expressed in the cell cycle and late cellular differentiation phase during adipogenesis. To investigate this dual role 
of pRb in the early and late stages of adipogenesis we used microarrays to perform a comprehensive systems-level 
analysis of the common transcriptional program of the classic 3T3-L1 preadipocyte cell line, wild-type mouse 
embryonic fibroblasts (MEFs), and retinoblastoma gene-deficient MEFs (Rb-/- MEFs). 

Findings: Comparative analysis of the expression profiles of 3T3-L1 cells and wild-type MEFs revealed genes 
involved specifically in early regulation of the adipocyte differentiation as well as secreted factors and signaling 
molecules regulating the later phase of differentiation. In an attempt to identify transcription factors regulating 
adipogenesis, bioinformatics analysis of the promoters of coordinately and highly expressed genes was performed. 
We were able to identify a number of high-confidence target genes for follow-up experimental studies. 
Additionally, combination of experimental data and computational analyses pinpointed a feedback-loop between 
Pparg and Foxol. 

To analyze the effects of the retinoblastoma protein at the transcriptional level we chose a perturbated system 
(Rb-/- MEFs) for comparison to the transcriptional program of wild-type MEFs. Gene ontology analysis of 64 
deregulated genes showed that the Rb-/- MEF model exhibits a brown(-like) adipocyte phenotype. Additionally, 
the analysis results indicate a different or additional role for pRb family member involvement in the lineage 
commitment. 

Conclusion: In this study a number of commonly modulated genes during adipogenesis in 3T3-L1 cells and MEFs, 
potential transcriptional regulation mechanisms, and differentially regulated targets during adipocyte differentiation 
of Rb-/- MEFs could be identified. These data and the analysis provide a starting point for further experimental 
studies to identify target genes for pharmacological intervention and ultimately remodeling of white adipose tissue 
into brown adipose tissue. 



* Correspondence: hubert.hackl@i-med.ac.at 

^Biocenter, Division of Bioinformatics, Innsbruck IVledical University, Austria 
Full list of author information is available at the end of the article 

O© 201 1 Hackl et al; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative ComrTions 
BiolVlGCl Centrsl Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 



Hakim-Weber et al. BMC Research Notes 201 1, 4:157 
http://www.biomedcentral.eom/1 756-0500/4/1 57 



Page 2 of 13 



Introduction 

Evidence accumulating during the past decades has 
convincingly revealed that adipocytes and adipose tissue 
are not only acting as a storage depot for fat, but are 
also actively involved in regulating whole body energy 
balance, and much attention has been dedicated to deci- 
pher molecular events and to identify factors involved in 
fat cell development. Although many molecular regula- 
tors have been identified, a comprehensive catalogue of 
components regulating adipogenesis is still missing. 
Therefore, a number of expression profiling studies 
using microarrays have been performed to monitor glo- 
bal gene expression profiles during in vitro adipocyte 
differentiation of several cell models [1]. 

Recently, we have studied global gene expression of 
3T3-L1 cells during adipogenesis and identifying mole- 
cular mechanisms, biological processes, key enzymes, 
including many characterized and uncharacterized tran- 
scriptional regulators, thereby providing a gene expres- 
sion atlas using functional annotation [2]. We showed 
that the mitotic clonal expansion phase plays an impor- 
tant role during adipocyte differentiation in the studied 
mouse cell model. One group of proteins - the pocket 
proteins including the retinoblastoma protein (pRb) and 
other members (pl30, pl07) - was expressed during 
both, the cell cycle phase and late cellular differentia- 
tion. This dual role of pRb in the early and late stages 
of adipogenesis was also described in other studies [3-6]. 
Pathways regulated by the retinoblastoma gene/protein 
might therefore be interesting novel targets in the com- 
bat to curb or to treat obesity and associated disorders. 

Several studies on mechanisms and interactions have 
elucidated how pRb - the product of the Rbl gene - fulfills 
its functions. pRb is phosphorylated by cyclin-dependent 
kinases thereby regulating E2F activity [5]. Other mechan- 
isms include co-factor activity for none-E2F transcription 
factors, recruitment of chromatin remodeling complexes 
(including HDACs) and action as a non-chromatin asso- 
ciated protein adaptor [4,6,7]. The first evidence that 
pocket proteins are involved in adipocyte differentiation 
emerged by showing interactions with the C/EBP tran- 
scription factors [8,9] . It was shown that the ability of the 
C-terminally truncated large T antigen (TAg) to inhibit 
differentiation of white preadipocyte cell lines is depen- 
dent on the ability to sequester pRb family members [10]. 
Further studies showed that pRb and pl07/pl30 may have 
opposite effects on adipocyte differentiation [11,12]. 

Rb-/- mouse embryonic fibroblasts (MEFs) are unable 
to undergo adipose conversion using standard adi- 
pogenic inducers [8,13,14]. It was suggested that the 
underlying mechanism did not directly mirror cell cycle- 
related aspects of pRb function, as Rb-/- cells enter and 
exit mitotic clonal expansion with the same kinetics as 
wild-type cells [13]. Further it was shown that pRb is 



required for preadipocyte differentiation in vivo [12,15] 
and functions as a molecular switch determining white 
versus brown(-like) adipocyte differentiation [14,16]. 
Of note, adipose tissue-specific inactivation of Rb 
prevents [15] and haploinsufficiency of Rb attenuates 
diet induced obesity and insulin resistance [17]. pRb 
might also be involved in the lineage commitment to 
osteogenic or adipogenic differentiation of mesenchymal 
stem cells [18]. 

Thus, given the accumulated evidence of the involve- 
ment of pRb in adipogenesis we initiated a large scale 
gene expression study to elucidate its role during adipo- 
genic differentiation. We first performed microarray 
gene expression analysis of wild-type mouse embryonic 
fibroblasts (MEFs) undergoing adipogenesis and com- 
pared these expression data with microarray gene 
expression data from 3T3-L1 cell adipogenesis to iden- 
tify a transcriptional program common between the two 
cell models. Additionally, bioinformatics analyses were 
carried on the promoters of highly expressed genes to 
identify potential transcription factor binding sites and 
to generate hypotheses on transcription factor binding. 
To analyze the effects of the retinoblastoma protein at 
the transcriptional level we chose a perturbated system 
(Rb-/- MEFs) for comparison to the transcriptional 
program of wild-type MEFs. 

Results and Discussion 

Common transcriptional program of adipocyte 
differentiation in MEFs and 3T3-L1 cells 

To compare transcriptional regulation during adipocyte 
differentiation of the established 3T3-L1 cell line and 
mouse embryonic fibroblasts (MEFs), gene expression in 
MEFs was studied in three independent experiments 
with whole genome cDNA microarrays (>27k elements) 
using seven time points (dO, dl, d2, d3, d4, d6, dlO) 
after hormonal induction in relation to expression levels 
at the preconfluent stage (similar to the generation of 
data from the previous 3T3-L1 study [2]) (Additional 
file 1). We validated the microarray data from MEFs by 
quantitative reverse transcription polymerase chain reac- 
tion (qPCR) analysis of the first differentiation experi- 
ment (see Figure SI in Additional file 2) with six genes 
at several time points. Applying linear regression of the 
resulting 18 data points showed a coefficient of determi- 
nation (7?^) of 0.92 (Pearson correlation coefficient 
r = 0.96), which represented good agreement between 
microarray and qPCR results. Replicate results were also 
in good concordance with log2-fold changes averaged 
over the three independent differentiation experiments. 
Differentially expressed genes with similar combined 
expression profiles of the 3T3-L1 and MEF cells were 
grouped into 6 clusters by k-means clustering (Figure 1 
and Additional file 1). Proficient quality of clustering 
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Figure 1 Gene expression levels (log2-fold change) of genes from different time points during adipocyte differentiation in 3T3-L1 and 
IVIEF cells. (A) Compact heatmaps from k-means clustering. (B) Line drawing centroid view (mean +/- SD). (C) Tliree-dimensional visualization 
of Principal Component Analysis (PCA) (color coding according k-means clustering) demonstrates cluster separation. (D) Line drawing for 
expression course of each gene individually. 
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was evident from different visualizations of the gene 
expression levels and separation of clusters using Princi- 
pal Component Analysis (PCA) in three-dimensional 
view of the first three principal components and color- 
coding according to k-means clustering (Figure 1). 
Functional assignment of genes within each cluster by 
gene ontology analysis revealed a number of noticeable 
biological processes regulated at the transcriptional level 
(Figure 2). Most significantly overrepresented biological 
processes were 'cell cycle', 'cell division', 'DNA replica- 
tion', and 'mitosis' indicating mainly the involvement of 
genes of cluster E and cluster B in cell cycle processes. 
Distinct peaks in the expression profiles of these genes 
at dl and d3 showed that MEFs underwent two cycles 
of mitotic clonal expansion comparable to the 3T3-L1 
cells showing - at the used sampling rate - one peak in 
the expression profiles at 24 h [2] which could reflect 
the first of two rounds of mitotic clonal expansions 
[19,20]. This observation is supported by reported phos- 
phorylation events of the pRb in 3T3-L1 cells during 
differentiation [3,13] and the cyclin A2 (Ccna2) expres- 
sion profile as shown in Additional file 1 and Additional 
file 2. 

As evident from the gene ontology analysis of genes 
in cluster E, significantly overrepresented biological 
processes were 'fatty acid metabolism', 'lipid metabo- 
lism', 'metabolism', and 'fatty acid beta-oxidation'. Well 
known genes important for the development of the 
adipocyte phenotype were combined in cluster F with 
increasing expression over time, including transcription 
factors and nuclear receptors Pparg, Cebpa, Srebfl, 
Xbpl, Thra, Statl, Nr0b2 (SHP), Nrldl (REV-ERB 
alpha), Nr3cl (GR), Foxcl, Foxol, Zhxl, and Zfp503. 
The recently identified gene coding for adipocyte 
plasma membrane-associated protein [21] was also 
shown to be highly expressed at late stages of adipo- 
genesis of MEFs (cluster F). 

Highly expressed genes over the differentiation course, 
especially at the middle stage of adipogenesis (after 
mitotic clonal expansion) could be found in cluster A. 
In this group we observed not only transcription factors 
like Klf9 [22], and Foxo3a [23] but also secreted factors 
(as studied in human [24]) including Apod [25], CcdcSO 
[26], Adamts5 [27], Mmp2 [28], and Ptn [29], as well as 
enzymes, signaling molecules, and other factors such as 
Xdh [30], Enpp2 [31], Rgs2 [32], Txnip [33], Tsc22d3 
[34], H6pd [35], Igfbp2 [36], and Igfbp3 [37]. Hence, the 
genes in cluster A represent potential novel players in 
adipogenesis and are therefore candidates to be studied 
in more detail experimentally. 

To identify candidates selectively associated with the 
MEF cell model, genes upregulated during MEF differ- 
entiation which showed marginal changes in 3T3-L1 
adipocytes compared to the preconfluent stage of 3T3- 



Ll are of major interest (cluster D). Genes in this clus- 
ter were overrepresented in the GO terms 'heart devel- 
opment', 'negative regulation of cell proliferation', 'cell 
adhesion', 'extra cellular matrix organization and biogen- 
esis', and 'positive regulation of angiogenesis'. The 
expressed transcription factors include Arnt, Arnt2, 
Atf5, Ets2, Egrl, Fosl2, JunB, Nfia, Prdm4, Stat2, Stat3, 
Tef, and Tcf20 which might be involved specifically in 
early regulation of adipocyte differentiation in the pri- 
mary cell model as opposed to the already committed 
3T3-L1 preadipocytes. 

Transcription factors regulating gene expression during 
adipocyte differentiation 

The terminal development of a fat cell is driven by a 
transcriptional cascade governed mainly by Pparg and 
Cebpa (as reviewed in [38,39]). The transcriptional 
events directing early adipogenic processes, concerted 
by molecular events prior to the activation of Pparg and 
Cebpa have not been exhaustively studied. Hence, we 
asked if there are common regulators of expression of 
genes within each cluster, especially those with expres- 
sion profiles similar to a transcription factor and its tar- 
get genes. To address this question, we calculated 
overrepresentation of sequence motifs (response ele- 
ments) within a genomic region (including promoter) 
from -3 kb to +2 kb relative to transcription start site 
(TSS) of those genes (Additional file 3). 

From the transcription factors with exclusively overre- 
presented binding sites (Figure 3A) in promoters of 
genes from cluster F (i.e. genes with increasing expres- 
sion towards terminal adipocyte differentiation) only 
Foxol and the key player in adipogenesis Pparg have 
expression profiles in cluster F. It has been previously 
shown that Foxol can interact with Pparg [40]. Further- 
more, Foxol can bind to the Pparg promoter region and 
suppress Pparg expression [41]. Our analyses suggest 
that a feed-back loop from Pparg to Foxol might exist. 
To provide further evidence we analyzed two datasets of 
DNA binding by Pparg identified by chromatin immu- 
noprecipitation (ChIP) followed by sequencing [42] or 
by microarray analysis [43] during 3T3-L1 adipocyte dif- 
ferentiation. Integrative analysis of these datasets and 
Pparg binding sites in genomic regions around genes 
with expression profile similar to that of Pparg in both 
the 3T3-L1 cells and MEFs (cluster F), revealed that 
Foxol might indeed be regulated by Pparg as indicated 
in Figure 3B. Thus, our data, experimental results from 
other studies, and further computational analysis pin- 
pointed a feed-back loop between Pparg and Foxol. 

We have recently identified 2310001A20Rik (Apmap), 
as well as the retrotransposed genes 6530401D17Rik 
(Arxesl) and 2900062LllRik (Arxes2) as Pparg targets 
and that they are required for adipogenesis [21,44]. 



Hakim-Weber et al. BMC Research Notes 201 1, 4:157 
http://www.biomedcentral.eom/1 756-0500/4/1 57 



Page 5 of 1 3 



0.05 


adj. p-value 


lE-20 




130 


'Ioel0(adj. 


2aQ 




3:0007049 
D:0051301 
3:0006260 
3:0007067 
3:0006631 
3:0006929 
3:0008152 
3:0006635 
3:0007059 
3:0000910 
3:0007507 
3:0008285 
3:0006928 
3:0000074 
3:0007001 
3:0008380 
; 0006397 
: 0rinF?-(57 
jO:i7I55 



J0ji:.334 
uuu6y86 
G': -0007017 
00:0006268 
00:0006323 
00:0016126 
00:0007050 
00:0007010 
00:0050885 
O'j 0006270 
■-rO 0008:83 



OU.OD00084 
50:0006099 
30:0030198 
50:0030833 
30:0001525 
30:0006281 
30:0043433 
30:0006118 
30:0006259 
00:0051258 
00:0009165 
00:0007266 
30:0046766 
30:0007015 
30:0008360 
30:0007031 
30:0006974 
30:0001503 
30:0009636 
30:0000059 
30:0044267 
30:0000122 
30:0006468 
30:0030036 
30:0006817 
30:0009058 
30:0006914 
30:0051017 
30:0008286 
30:0007411 
30:0006096 
30:0046892 
30:0043066 
30:0008219 
30:0008217 
30:0006694 
30:0008610 
00:0051028 
30:0007265 
00:0008652 
00:0009408 
00:0016358 
00:0007018 
Oi:i.0008814 
00 0007259 
O'J 0000074 

03 

-11 

Ou:UUU6yl5 
30:0030087 
30:0007417 
30:0007626 
30:0006979 
30:0007160 
30:0000398 
30:0008015 
30:0001701 
30:0000160 
30:0042254 
30:0016477 
30:0001764 
30:0006836 
30:0006139 
30:0006520 
O'j 0045786 
>'"' j0j8:84 



Oij .0043087 
e 0:0006350 
G 0:0042981 



cell cycle 
cell division 
DNA replication 
mitosis 

fatty acid metabolism 
lipid metabolism 
metabolism 

fatty acid beta-oxidation 

chromosome segregation 

cytokinesis 

heart development 

negative regulation of cell proliferation 
cell motility 

regulation of progression through cell cycle 

chromosome organization and biogenesis(sensu Eukaryota) 

RNA splicing 

mRNA processing 

protein folding 

cell adhesion 

nucleosome assembly 

response to unfolded protein 

microtubule-based process 

DNA unwinding during replication 

DNA packaging 

sterol biosynthesis 

cell cycle arrest 

cytoskeleton organization and biogenesis 

regulation of balance 

DNA replication initiation 

cell proliferation 

regulation of cell growth 

rRNA processing 

S phase of mitotic cell cycle 

tricarboxylic acid cycle 

extracellular matrix organization and biogenesis 
regulation of actin filament polymerization 
angiogenesis 
DNA repair 

negative regulation of transcription factor activity 

electron transport 

DNA metabolism 

protein polymerization 

nucleotide biosynthesis 

Rho protein signal transduction 

positive regulation of angiogenesis 

actin filament organization 

regulation of cell shape 

peroxisome organization and biogenesis 

response to DNA damage stimulus 

ossification 

response to toxin 

protein import into nucleus, docking 
cellular protein metabolism 

negative regulation of transcription from RNA polymerase II promoter 

protein amino acid phosphorylation 

actin cytoskeleton organization and biogenesis 

phosphate transport 

biosynthesis 

autophagy 

actin filament bundle formation 
insulin receptor signaling pathway 
axon guidance 
glycolysis 

negative regulation of transcription, DNA-dependent 
negative regulation of apoptosis 
cell death 

blood pressure regulation 

steroid biosynthesis 

lipid biosynthesis 

mRNA transport 

Ras protein signal transduction 

amino acid biosynthesis 

response to heat 

dendrite development 

microtubule-based movement 

sodium ion transport 

JAK-STAT cascade 

regulation of progression through cell cycle 
RNA processing 
cholesterol metabolism 
cellular process 

generation of precursor metabolites and energy 
anterior/posterior pattern formation 
apoptosis 
hemopoiesis 

central nervous system development 

locomotory behavior 

response to oxidative stress 

celUmatrix adhesion 

nuclear mRNA splicing, via spliceosome 

circulation 

embryonic development (sensu Mammalia) 
two-component signal transduction system (phosphorelay) 
ribosome biogenesis and assembly 
cell migration 
neuron migration 
neurotransmitter transport 

nucleobase, nucleoside, nucleotide and nucleic acid metabolism 
amino acid metabolism 

negative regulation of progression through cell cycle 

positive regulation of cell proliferation 

cell-cell adhesion 

regulation of OTPase activity 

transcription 

regulation of apoptosis 



Figure 2 Gene ontology analysis. Heatmap for overrepresented gene ontology terms for genes in 6 clusters over several time points during 
adipocyte differentiation in 3T3-L1 and MEF cells as visualized in Figure 1 (color coding according to legend at the top, only gene ontology 
terms with adjusted p < 0.05 were considered) is shown. 
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Figure 3 In silico promoter analysis and genomic organization. (A) Statistically overrepresented (p < 0.05) potential transcription factor 
binding sites (TFBS) identified by in silico analysis of genomic regions from -3 kb to + 2 kb around transcription start sites of genes with 
increasing gene expression at late stages of adipocyte differentiation (cluster F) using position frequency (weight) matrices from TRANSFAC and 
JASPAR. Corresponding sequence logos are shown and number of targets within the dataset (cluster) is given and tested against the number of 
targets from all RefSeq transcripts using Fishers exact test. Note only transcription factor motifs are considered which are exclusively 
overrepresented in cluster F. (B) The genomic organization around Foxol using the UCSC genome browser (NCBI37/mm9) with high similarity 
to the Pparg-Rxra DR-1 motif (similarity score > 0.90) (blue strand, red '-'strand), Pparg binding regions as identified by previous studies at day 
6 [42] (dark red) and day 10 [43] (orange) of 3T3-L1 adipocyte differentiation. 



A number of other genes annotated as Riken clones (like were Acadm, Adcy5, Dgatl, Etfb, Echl, Hmoxl, Mdhl, 

5730469M10Rik and F630110N24Rik) included in cluster MgstS, Pparg, Pla2gl5, Phldbl, Sorbsl, and Srebfl as 

F show noticeable binding by Pparg in 3T3-L1 cells in indicated in both chromatin immunoprecipitation 

the first intron and hence might be interesting targets in (ChIP) studies [42,43] of Pparg binding at the same geno- 

follow-up experimental studies. Additional target genes mic region near TSS (< 1 kb) or in an intron (peak 
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height > 70) (see additional website http://icbi.at/adipo). 
There was also supporting evidence of overrepresented 
known cell cycle related transcription factor binding sites 
in the genes of cluster B, C, and E (including NF-Y, E2F, 
cMyc:Max, USF, c-ETS). In cluster A exclusively overre- 
presented transcription factor binding sites includes 
STATS, STAT5A, STAT5B, MEF-2, Handl:E47, and NF- 

1. The genes encoding for STATS and NF-I could also 
be found in this cluster, indicating that these two tran- 
scription factors share similar expression profiles as their 
target genes in cluster A. 

Transcriptional program of Rb-/- MEFs 

Recent studies showed that the retinoblastoma protein 
(pRb) was expressed in the cell cycle and late cellular 
differentiation phase during adipogenesis. The interest- 
ing Rbl gene expression profiles during ST3-L1 and 
wild-type MEF adipogenesis are also shown in Addi- 
tional filel. Consequently, we studied the effect of per- 
turbation of the MEF adipogenesis system by lack of 
Rbl by a direct comparison of the global expression 
levels of Rb-/- MEF (MES) cells to those of the wild- 
type MEFs (MEFA) at 3 time points (dl, dS, d8) during 
adipocyte differentiation in relation to the expression at 
day 0, immediately before hormonal induction using 
again whole genome cDNA microarrays (> 27k ele- 
ments). Since Rb-/- MEFs do not differentiate into adi- 
pocytes in response to a standard adipogenic cocktail, 
the proadipogenic Pparg agonist Rosiglitazone was 
added in both the Rb-/- MEF (MES) and the Rb+/ + 
MEF (MEFA) differentiation experiments. We calculated 
several distance measures (sum of differences, Manhat- 
tan distance, Euclidean distance, and Pearson correla- 
tion) between the expression profiles of MES and MEFA 
cells for S118 differentially expressed ESTs (Additional 
file 4) and focused on 64 upregulated genes (sum of dif- 
ferences > 4) and 90 downregulated genes (sum of dif- 
ferences < -4). For clarity, gene expression levels of only 
the most deregulated ESTs were visualized as heatmap 
in Figure 4 (2S ESTs with sum of differences > 5 and 41 
ESTs with sum of differences < -5). 
As evident from previous studies, microarray and qPCR 
analysis revealed that Ucpl, Ppargcla, Prdml6, and 
Cidea are highly upregulated in Rb-/- MEF (MES) adi- 
pocyte differentiation compared to Rb+/+ MEF (MEFA) 
adipocyte differentiation (see [14,16,45], Additional file 

2, and Additional file 4). Gene ontology analysis from 
the 64 deregulated genes corroborated the brown(-like) 
adipocyte phenotype of the Rb-/- MEF model: the most 
significant overrepresented gene ontology term was 
'generation of precursor metabolites and energy' (FDR = 
3.1E-05) (including electron transferring flavoprotein, 
alpha polypeptide (Etfa) and NADH dehydrogenase (ubi- 
quinone) Fe-S protein 1 (Ndufsl) also involved in 



oxidative phosphorylation a prominent process in brown 
adipose tissue). Analysis with the tool ClueGO (Figure 
5) showed prominent gene ontology terms 'mitochon- 
drion organization', 'electron transport', 'triglyceride 
metabolic processes', and 'apoptosis'. These are related 
to hallmarks of a brown adipocyte phenotype such as 
mitochondrial biogenesis, metabolic processes, electron 
transport, oxidative phosphorylation, and uncoupling. 
Many of the activated genes are encoding for enzymes 
involved in glucose or lipid metabolism and one gene 
for a transcription factor (Myc). 

It is noteworthy that gene set enrichment analyses of 
the ranked gene list with gene sets from brown adipose 
tissue versus white adipose tissue [46] and brown and 
white preadipocytes at the differentiating stages [47] 
showed a high enrichment score at least to the first 
gene set supporting the notion that Rb is involved in 
the switch between the white and brown adipocyte dif- 
ferentiation process but also indicated expression phe- 
notypic differences between brown adipocytes and 
differentiated Rb-/- MEFs (Additional file 5). These 
results suggest a different and/or additional role of pRb. 
We have seen that Rb-/- MEFs behave more like pri- 
mary cells from inguinal white adipose tissue than clas- 
sic brown adipose tissue [48]. It was also shown that 
pRb could have a different role as it might be involved 
in the lineage commitment to either osteogenic or adi- 
pogenic differentiation of mesenchymal stem cells [18]. 

Conclusion 

In this work we studied adipocyte differentiation of 
wild-type MEFs and Rb-/- MEFs into mature adipocytes 
using microarrays and comprehensive bioinformatics 
analyses. We identified novel candidates involved in 
early adipogenesis as well as high-confidence genes tar- 
geted by Pparg. Our findings corroborate the brown-like 
phenotype of the Rb-/- MEFs. Based on the analyses we 
also propose the existence of a feed-back loop between 
Pparg and Foxol adding a novel facet to the regulatory 
network orchestrating adipocyte differentiation. 

The data and the analysis provide a starting point for 
further experimental studies to identify target genes for 
pharmacological intervention and ultimately remodeling 
of white adipose tissue into brown adipose tissue. 

Materials and methods 

Adipocyte differentiation of wild-type (Rb+/+) and Rb-/- 
mouse embryonic fibroblasts 

Generation of Rb-/-MEFs and wild-type MEFs were pre- 
viously performed at the Division of Cancer Biology, 
Danish Cancer Society, Copenhagen, DK and described 
in [49]; detailed information on the origin of these cells 
and disruption of the Rb-1 gene in mice can be found 
in [50]. A batch of the Rb-/- and Rb+/+ MEFs were 
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solute earner family 16 (monocarboxylic acid transporters), member 1 (Slc16a1) 
cell death-inducing DNA fragmentation factor, alpha subunit-like effector A (Cidea) 
UTP18, small subunit(SSU) processome component, homolog (yeast) (Utp18) 
ATG13 autophagy related 13 homolog (S. cerevisiae)(Atg13) 
0 protein-coupled receptor 176 (Opr176) 
lumican (Lum) 

superoxide dismutase 2, mitochondrial (Sod2) 
embryonic ectoderm development (Eed) 
BCL2/adenovirus E1B interacting protein 1. NIP2(Bnip2) 
aldolase 3, C isoform (Aldoc) 
nephronectin (Npnt) 
glypican 1 (Gpc1) 

glutamic pyruvate transaminase (alanine aminotransferase) 2 (Gpt2) 
cytoplasmic polyadenylation element binding protein 3 (Cpeb3) 
cellular retinoic acid binding protein II (Crabp2) 
metallothionein 1 (Uti) 
CD24a antigen (Cd24a) 
DNA-damage inducible transcript 3 (Ddit3) 

solute carrier family 1 (glial high affinity glutamate transporter), members (Slc1a3) 
pellino 1 (Peli1) 

RIKEN cDNA 1ig0002N15 gene (1190002N1SRik) 
glycogenin (Oyg) 

phosphodiesterase 3B. cGMP-inhiblted (PdeSb) 

fibroblast growth factor 10 (FgflO) 
tight Junction protein 1 (Tjpl) 

protein phosphatase 2C. magnesium dependent, catalytic subunit(Ppm2c} 
expressed sequence AW124722 (AW124722) 
regulator of G-protein signaling 12 (Rgs12) 
R-spondin 2 homolog (Xenopus laevis) (Rspo2) 
activated leukocyte cell adhesion molecule (Alcam) 
insulin- 1 ike growth factor 2 (Igf2) 
collagen, type I. alpha 2 (Col1a2) 
syndecan 2 (Sdc2) 
contactin 4(Cntn4) 
cell adhesion molecule 1 (Cadml) 
transmembrane 4 superfamily member 1 (Tm4$f1) 
spire homolog 1 (Drosophila)(Spire1) 
thioredoxin interacting protein (Txnip) 
pogo transposable element with KRAB domain (Pogk) 
sulfatase 2 (Sulf2) 

SERTA domain containing 4(Sertad4) 
chemokine (C-X3-C motif) ligand 1 (Cx3cl1) 
PERP. TP53 apoptosis effector (Perp) 
PQ loop repeat containing (Pqlc3) 
HI histone family, member 0 (HlfD) 
nuclear mitotic apparatus protein 1 (Numal) 
zinc finger. DHHC domain containing 8 (ZdhhcS) 
colony stimulating factor 1 (macrophage) (Csfl) 
exostoses (multiple)- 1 ike 3 (Extl3) 
p ro g ressive a nkyi osis (Ank) 

activated leukocyte cell adhesion molecule (Alcam) 
fibulin 2 (Fbln2) 

collagen, type XVIII. alpha 1 (Col18a1) 
sulfatase 2 (Sulf2) 

serine (or cysteine) peptidase inhibitor, clade B. member la (Serpinbia) 
RIKEN cDNA 4033427014 gene (4e33427D14Rik) 
four and a half LIM domains2 (Fhl2) 
collagen, type XVIII. alpha 1 (CollSal) 
angiotensin II receptor, type 2 (Agtr2) 
pleiotrophin (Ptn) 
pleiotrophin (Ptn) 
PR domain containing 4(Prdm4) 
matrix Gla protein (Mgp) 
talin 1 (TInl) 

Figure 4 Gene expression profiles in Rb-/- MEFs. Gene expression levels (log2-fold change) shown as heatmap for samples from 3 different 
time points (d1, d3, d8) related to dO immediately before hormonal induction of Rb-/- MEF (ME3) and Rb+/+ MEF (MEFA) differentiation. Only 
genes with sum of differences > 5 or sum of differences < -5 are shown. 




kindly provided by Jiri Lukas as already used in other 
studies (e.g. [14]). No animals or human subjects were 
directly involved in the current study that would make 
ethical approval and consent required. Wild-type mouse 
embryo fibroblasts (MEFs) were differentiated at the 
Department of Biochemistry and Molecular Biology, 



University of Southern Denmark, Odense, DK according 
to the protocol described in [13]. Briefly, cells were 
induced to differentiation by the MDI protocol (1 |iM 
dexamethasone (DEX), 0.5 mM isobutyl methylxanthine 
(MIX), and 5 \xg/m\ insulin). DEX and MIX were 
omitted after day 2, but insulin was added throughout 
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Figure 5 ClueGO analysis [57]of 64 genes upregulated during Rb-/- MEF adipocyte differentiation in comparison to Rb-i-/-i- MEF 
adipocyte differentiation (sum of differences of Iog2-fold change > 4). Nodes (circles) represent gene ontology terms. Connections 
between two nodes (edges) indicate that two gene ontology terms share genes from the considered dataset (agreement measure kappa > 0.3). 
The calculated kappa score is also used for defining functional groups, which are indicated by the same color. The most prominent gene 
ontology term for each group (theme) is highlighted: mitochondrion organization (red), eiectron transport (blue), trigiyceride metabolic process 
(green), apoptosis (orange). 



the differentiation. Cells were passaged and differen- 
tiated in AmnioMax (Gibco) supplemented with 7,5% 
fetal bovine serum. Nutrition media were changed every 
second day. Rb-/- MEFs and wild-type MEFs (Rb+/+) 
for comparison were differentiated using the identical 
protocol, but adding Pparg ligand BRL49653 (0.5 \iM) to 
the induction cocktail throughout the whole differentia- 
tion process. All differentiation experiments were per- 
formed in triplicate. Cells were harvested and total RNA 
was isolated at the preconfluent stage and at seven time 
points after induction (0, 1, 2, 3, 4, 6, and 10 days) with 
TRIzol reagent (Invitrogen). For the Rb-/- and Rh+/+ 
MEF differentiation, cells were harvested immediately 
before (reference) and at three time points (1, 3, and 8 
days) after induction. For each independent experiment 
RNA was pooled from three different culture dishes for 
each time point and from twenty dishes for each 



reference. The quality of the RNA was checked using 
Agilent 2100 Bioanalyzer RNA assays (Agilent Technol- 
ogies) by inspection of the 28S and IBS ribosomal RNA 
intensity peaks. 

Gene expression and data analysis 

The mouse cDNA microarrays and the hybridization 
protocols used have been described previously [2]. 
Briefly, the spotted microarrays contain 27,000 elements 
with mouse cDNA clones representing 16,000 different 
genes (UniGene clusters). 10-25 micrograms (dependent 
on replicated experiment) of total RNA from each time 
point was reverse transcribed into cDNA, which was 
then indirectly labeled with Cy5 and Cy3, respectively. 
To account for technical variation, procedures were 
repeated using the same samples with reversed dye 
assignment. The microarrays were pre-hybridized with 
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5x SSC, 0.1% SDS,1% BSA. Pair-wise labeled cDNA 
samples were combined and 20 \ig of mouse Cotl DNA 
and 20 i^g of poly(A) DNA were added. The mixture 
was hybridized onto the slides overnight at 42°C. Fol- 
lowing washing, slides were scanned with a GenePix 
4000B microarray scanner (Axon Instruments) at 10 i^m 
resolution. The resulting TIFF images were analyzed 
with GenePix Pro 4.1 software (Axon Instruments). Fea- 
tures were filtered for low-quality spots. Following sub- 
traction of the local background, the arrays were global 
median and dye-swap normalized using ArrayNorm [51] 
and the resulting ratios log2 transformed. All experi- 
mental parameters, raw data, and transformed data were 
uploaded to the microarray database MARS [52] and 
submitted via MAGE-ML export to a public repository 
(ArrayExpress [53], Accession Nos. E-MARS-12/E- 
MARS-13 and A-MARS-7). After outliers were removed 
the median expression values of replicated ESTs on the 
microarray were calculated. Only ESTs with z > 1.5 and 
significant replication {p < 0.05) in at least one time 
point using the z test and implying standard normal dis- 
tribution were considered differentially expressed and 
used for further analysis [54,55]. Expression values and z 
scores for differentially expressed ESTs were averaged 
over biological replicates and annotated using the 
RefSeq database (MegaBLAST E-value < lE-30) and 
Entrez Gene database. Expression profiles for MEF data 
and 3T3-L1 data [2] were combined based on the same 
EST ID and profiles were grouped into 6 clusters by k- 
means clustering and validated by principal component 
analysis. Gene ontology term overrepresentation analysis 
for the different gene sets (clusters) were performed 
using Fisher exact test and p-values were adjusted for 
multiple hypothesis tests based on the Benjamini-Hoch- 
berg method (http://genome.tugraz.at/ORA). Several dis- 
tance measures between the expression profiles of the 
Rb-/- MEF (ME3) and Rh+/+ MEF (MEFA) differentia- 
tion models (sum of differences, Manhattan distance, 
Euclidean distance, Pearson correlation coefficient) were 
calculated. Over-represented gene ontology (GO) terms 
for up- (sum of differences > 4) or downregulated (sum 
of differences < -4) genes were identified using DAVID 
[56] based on gene symbols. Up-regulated genes were 
functionally grouped into gene ontology networks using 
the Cytoscape plug-in ClueGO [57]. All calculations 
were implemented in PERL 5.8.0 or R 2.10.1, cluster 
analyses and visualizations were performed using Gen- 
esis [58]. Gene set enrichment analysis (GSEA) [59] was 
performed using a tool from http://www.broadinstitute. 
org/gsea. 

qPCR analysis 

Microarray expression results were confirmed with 
qPCR. cDNA was synthesized from 1.25 ^g total RNA 



in 20 |il using random hexamers and Superscript II 
reverse transcriptase (Invitrogen). The design of LUX^"^ 
primers was done using the Invitrogen web service. 
Quantitative RT-PCR analyses were performed starting 
with 12 ng reverse transcribed total RNA, with 0.5 x Pla- 
tinum Quantitative PGR SuperMix-UDG (Invitrogen), 
with a ROX reference dye, and with a 200 nM concen- 
tration of both LUX™ labeled sense and antisense pri- 
mers (Invitrogen) in a 25 \il reaction. Ribosomal 18S 
RNA amplifications were used to account for variability 
in the initial quantities of cDNA. The relative quantifi- 
cation with respect to the calibrator (preconfluent stage) 
was determined using the AACt method and compared 
with the normalized ratios resulting from microarray 
experiments. To better characterize Rb-/- MEF (ME3) 
and Rb+/-i- MEF (MEFA) cells, triplicate-qPCR analyses 
were performed with 4.5 ng cDNA, 200 nM primer, and 
Platinum SYBR green mastermix (Invitrogen). Data were 
analyzed using a real-time PGR management and analy- 
sis system [60]. Values were normalized to Uxt, Ct 
values and PGR efficiencies were derived as in [61]. 
All qPCR analyses were performed on an ABI PRISM 
7000 detection system (Applied Biosystems). Primer 
sequences can be found in Additional file 2. 

In silico analysis of transcription factor binding sites 

Statistically overrepresented (p < 0.05) potential tran- 
scription factor binding sites (TFBS) were identified by 
in silico analysis of genomic regions from -3 kb to -t- 2 
kb around transcription start sites of genes within each 
cluster and all RefSeq transcripts using known (TRANS- 
FAG [62], JASPAR [63]), and newly compiled position 
frequency (weight) matrices (http://genome.tugraz.at/ 
Logo). Potential binding sites were identified if the simi- 
larity score for each PWM was equal or above the 
threshold (allowing 1 binding site with lOkb of mouse 
coding sequences) based on an implementation of the 
Matlnspector algorithm [64]. The number of targets 
(showing at least one TFBS) within the gene set (cluster) 
was tested against the number of targets from all RefSeq 
transcripts using Fishers exact test and p-value was 
adjusted for multiple hypothesis testing based on the 
false discovery rate according to the method of Benja- 
mini-Hochberg [65]. Genome organization around genes 
(genomic region from 20 kb upstream of transcription 
start to 20 kb downstream of transcription end) showing 
similar expression profile as Pparg were constructed 
using UCSC genome browser (NCBI37/mm9), custom 
track data from chromatin immune precipitation 
followed by sequencing at day 6 [42] or microarray ana- 
lysis at day 10 [43] during 3T3-L1 adipocyte differentia- 
tion, and Pparg-Rxra direct repeats 1 (DRl) motifs 
(similarity score > 0.90) (see additional website http:// 
icbi.at/adipo). 
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